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Abstract 


The  differential  equation  u"(k,x)  +  [k  -V(xj]  u(k,x)  =  0  is  considered 

in  the  interval  -oo  <  x  <  oo.  From  a  knowledge  of  the  reflection  coefficient 

ikx 
b(k)  associated  vrith  a  wave  e    incident  from  the  left  it  is  shown  how  one 

can  calculate  the  function  V(x)  of  the  differential  equation  (when  V(x)  -  0 

for  X  less  than  some  constant).  The  method  used  is  due  to  I.M,  Gelfand  and 

B.M.  Levitan . 

Sufficient  conditions  such  that  a  function  b(k)  can  be  the  reflection 

coefficient  corresponding  to  a  differential  equation  of  this  form  are  given, 

2   -1 
An  example,  in  which  b(k)  »  -(k  +1)  ,  is  worked  out  in  detail.   It 

is  also  shown  how  one  can  obtain  V(x)  explicitly  for  a  large  class  of  reflec- 
tion coefficients. 

The  differential  equation  u"(k,x)  +  k^n^(x)  u(k,x)  =  0  in  the 
interval  -oo  <  x  <  oo  is  also  considered,  the  problem  being  to  determine  the 
function  n(x)  from  a  knowledge  of  the  reflection  coefficient,  assuming  that 
n(x)  ■  constant  for  x  less  than  some  given  value. 
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1,  Introduction 

It  is  a  customary  procedure  in  mathematical  physics  to  make  predic- 
tions about  particular  physical  situations  on  the  basis  of  a  general  theory 
and  then  to  compare  those  predictions  with  experimental  data.  This  procedure 
is,  of  course,  necessary  to  check  the  validity  of  a  new  theory.  But  tihen   a 
physical  theory  has  evoked  sufficient  confidence  because  of  its  continviing 
success  over  many  years  one  need  no  longer  look  upon  such  a  procedure  as  a 
test  of  the  theory.  Rather,  the  theory  becomes  a  tool  by  means  of  which 
one  can  describe  in  detail  the  various  special  cases  arising  when  particular 
conditions  are  added  to  the  general  theoretical  assumptions.  Experimental 
data  for  a  special  case  are  then  useful  to  measure  or  make  determinate  those 
physical  quantities  of  the  general  theory  ^ich  otherwise  are  arbitrary.  In 
the  application  of  a  trusted  theory  to  special  cases  it  still  seems  to  be 
customary,  however,  to  proceed  in  the  same  way  as  in  the  verification  of  a 
new  theory.  One  makes  certain  assumptions  about  the  relevant  physical  sit- 
uation, one  predicts  the  consequences  of  these  assumptions  on  the  basis  of 
the  general  theoiTr,  and  finally  one  makes  an  eB?)irical  check  of  the  assun^)- 
tions  by  means  of  experimental  data.  The  assxuaptions  can  be  corrected  accord- 
ing to  the  experimental  results  and  the  procedure  repeated  until  the  desired 
accuracy  is  reached.  In  the  case  of  an  old  and  tested  theory,  however,  it  would 
appear  more  logical  to  adopt  a  different  procedure,  if  possible,  A  gain  in 
time,  accuracy  and  understanding  might  be  expected  if  one  were  able  to  recon- 
stnict  directly  from  the  experimental  data  the  physical  conditions  correspond- 
ing to  a  special  case  of  the  general  theory. 

The  problem  of  making  srch  a  reconstruction  is  the  direct  inverse  of 
the  problem  ordinarily  undertaken.  This  reconstruction  is  closely  related  to 
the  problem  of  synthesis  wherein  the  task  is  to  construct  special  physical  con- 
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ditions  which,  according  to  the  general  theory,  will  produce  experimental 
data  given  beforehand. 

An  example  of  the  'inverse  problem'  is  the  attempt  to  determine 
the  ionization  density  of  the  ionosphere  from  the  time  delay  of  a  pulse 
radio  wave  has  been  transmitted  from  the  earth  and  relected  back  to  the 
earth  by  the  ionosphere.  The  synthesis  problem  is  a  common  one  for  elec- 
tical  engineers  who  are  often  required  to  design  networks  which  produce  a 
prescribed  response  to  a  prescribed  input  voltage  over  a  given  frequency  band. 

Some  work  has  been  done  in  recent  years  on  inverse  problems  suggest- 
ed by  quantum  mechanics  and  associated  with  the  differential  equation 

-^  u(k,x)  +  fk^-  V(x)l  u(k,x)  -  0  . 

For  the  quantum  mechanical  application  the  boundary  condition  u(k,0)  ■  0  is 
assiuned  and  the  variable  x  lies  in  the  interval  0  <  x  <   oo.     The  given  experi- 
mental data  provide  a  kr:owled>j,e  of  the  phase  fiu-iction  0(k)  vrtiich  appears  in 
the  asymptotic  form  of  u(k,x)  for  large  positive  x: 


u(k,x)  r>^  A(k)   sin  [kx  +  0(k)], 


+    00 


The  problem  is:  given  0{k)}   determine  V(x),  The  uniqueness  of  the  solution 

^2 
when  the  operator  -  — r  +  V(x)  has  no  bound  states  (proper  eigenf unctions) 

dx*^ 

was  proved  by  N.  Levinson'-  -^  and  V.A.  Marcenko'-  -' .  A  similar  inverse  problem 
for  the  same  differential  equation  arises  when  a  finite  interval  is  chosen 


and  homogeneous  boundary  conditions  are  applied  at  both  ends:  given  the  two 

^2 
sets  of  eigenvalues  of  the  operator  -  — ^  •»•  V(x)  corresponding  to  two  different 

dx^ 
pairs  of  boundary  conditions;  determine  V(x),  The  uniqueness  of  V(x)  in  this 

situation  was  proved  by  N,  Levinson'-^-'  and  G.  Borg'-  -^  and  V,A.  Marcenko'-  -' » 

I.M.  Gelfand  and  B.M.  Levitan"-  ^  have  shown  how  to  construct  the  quantity  V(x) 
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for  the  interval  0  <  x  <  oo  from  the  spectral  measure  function  associated  with 

,2  " 
the  operator  -  — r-  +  V(x)  and  boundary  conditions.  Another  solution  by  an  en- 

dx"^ 
tirely  different  method,  capable  of  including  the  case  in  which  bound  states 

exist,  was  given  by  W.  Kohn  and  R.  Yost L  J,  M,Krein'"j  has  solved  the  inverse 

problem  in  a  different  way  for  the  finite  interval, 

A  very  clear  English  siimmary  of  the  paper  by  I.M,  Gelfand  and  B.M. 
Levitan'^-',  with  an  additional  note  on  the  error  due  to  a  small  perturbation 
in  the  spectral  measure  function,  has  been  published  by  N.  Levinson'- -' . 

It  will  be  the  purpose  of  this  paper  to  treat  the  same  ordinary  differ- 
ential equation  in  the  interval  -oo  <  x  <  oo.  The  information  assumed  given 
will  be  the  reflectLoi  coefficient  b(k)  when  the  solution  u(k,x)  of  the  differ- 
ential equation  has  the  asymptotic  forms 

u(k,x)  ">   e^^  +  b(k)  e"^^  ,         x  -»•  -oo, 
u(k,x)  -»>  a(k)  e   ,  x  -»•  +oo  . 

The  function  V(x)  will  be  obtained  from  b(k)  by  an  adaptation  of  the  method 
of  I.M.  Gelfand  and  B.M.  Levitan '--'-J ,  and  an  explicit  example  will  be  worked 
out.  Moreover,  sufficient  conditions  will  be  derived  for  a  function  b(k)  to 
be  the  reflection  coefficient  of  such  a  differential  equation  (this  is  necess- 
ary for  the  synthesis  problem),  and  a  method  for  solving  a  large  class  of  in- 
verse problems  explicitly  will  be  given.  Finally,  the  differential  equation 

2 
-^  u(k,x)  +  k^  n^(x)  u(k,x)  «  0, 
dx*^ 

for  which  the  problem  is  to  find  the  function  n(x),  will  be  considered. 

The  solutions  obtained  in  this  article  appear  to  have  n\unber  of  appli- 
cations. The  synthesis  problems  for  a  non-uniform  transmission  line  (in  which 
either  the  inductance  or  the  capacitance  is  constant)  and,  under  certain  con- 
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ditions,  for  a  waveguide  of  variable  cross-section  can  be  treated  by  the 
method  described  here.  A  synthesis  problem  in  which  a  variable  stratified 
dielectric  is  to  be  constructed  so  that  the  reflection  coefficient  of  an  electro- 
magnetic wave  will  vary  in  a  prescribed  way  with  frequency  can  also  be  treated. 
Another  synthesis  problem  involves  a  plane  wave  incident  on  a  plane  interface 
of  a  variable  dielectric  slabj  the  reflected  wave  amplitude  variation  as  a  func- 
tion of  the  angle  of  incidence  is  prescribed,  and  the  variation  of  the  dielectric 
properties  of  the  slab  can  be  found  so  as  to  produce  the  prescribed  reflected 
wave  amplitude  variation.  The  problem,  already  mentioned,  of  determining  the 
ionization  density  of  a  stratified  ionosphere  is  another  possible  application. 
In  this  investigation  the  preliminary  information  used  concerning  the 
differential  equation 


dx2 


u(k,x)  +  fk^-  V(x)  u(k,x)|  -  0 


and  the  general  facts  about  the  spectral  theory  and  about  the  scattering  matrix 
associated  with  the  differential  equation  in  the  interval  -oo  <  x  <  oo  can  be 
fo\md  in  the  Appendix  17  of  K.  Friedrich's  notes  on  the  spectral  theory  of  linear 
operators  in  Hilbert  space"-  -■ . 
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2,  The  basic  differential  equation  and  some  of  Its  properties 

If  a  continuously  differentiable  wave  function  u(k,x)  satisfies  the 
differential  equation 


(2.1)  i^i^  * 


k^-  V(x) 


u(k,x)  -  0, 


vrtiere  the  function  V(x)  is  real,  bounded,  piecewise  continuous  and  belongs 
to  L,  in  the  interval  -oo  <  x  <  oo,  then  it  is  well  known  that  the  asymptotic 
behavior  of  u(k,x)  can  be  expressed  in  the  form 

u(k,x)  ^  Y^(k)  e^^+  5^(k)  6*^'°',     X  -9.  +00 
u(k,x)  ^   Y.(k)  6^*™+  5_(k)  e    ,     x  -»►  -oo  . 

The  quantities  y    ,   5       are,  as  indicated,  functions  of  k  alone  and  independ- 
ent of  X.  Moreover,  a  linear  relation  exists  between  the  components  of  the 


expressed  in  terms  of  a  two  by  two  square  matrix,  the  so-called  scattering  matrix 

St  I  ^  I   "Sic  /  .  The  four  components  of  S  are  functions  of  k  alone  which 
are  the  saiM  for  all  solutions  u(k,x)  of  (2.1)  and  are  fixed  by  the  function  V(x) 
appearing  in  (2.1).  If  we  write  S  in  the  form 

(a(k)  a'(k)\ 
p(k)  p  (k)  / 

the  following  relations  between  the  components  of  S  hold  for  real  k: 

o(-k)  -  aTU  -  a'(k) 
(2.U)  P(-k)  -  WT'   p'(k) 
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From  the  last  equation  in  (2.U)  it  follows  that 

(2.5)  la|2>  |p|2;   |a|2>i, 

for  finite  values  of  |a|  and  |p|  when  k  is  real.  It  is  also  known  that  a(k)  and 
p(k)  are  analytic  functions  of  k  in  the  complex  plane. 

It  is  possible  to  find  a  solution  of  (2.1)  which  has  the  following  asymp- 
totic behavior: 

ti(k,r)  -  e^'"  +  b(k)e"^^,     X  ->►  -00 

-  a(k)  e         ,      X  -*  +00. 

From  the  definition  of  the  scattering  matrix  we  have 

(2.7)  b(k)  -  |{|}  ,  a(k)  -  ^    , 

and  from  (2.U)  it  follows  that 

(2.8)  |a|2  +  |b|2  -  1 

when  k  is  real.  In  accordance  with  the  asymptotic  behavior  of  u(k,x)  given  by 

(2.6)  the  quantity  b(k)  is  called  the  reflection  coefficient  and  the  quantity 

a(k)  the  transmission  coefficient, 

^2 
In  general,  for  the  interval  -co  <  x  <  +00,  the  operator  -  — r-  ••■  V(x) 

dx 
of  equation  (2.1)  has  a  spectral  representation  with  a  spectrum  composed  of  a 

discrete  part  and  a  continuous  part.  The  discrete  spectrum  has  for  its  corres- 
ponding eigenf unctions  those  solutions  of  (2.1)  which  belong  to  Lp.  For  values 
of  k  in  the  discrete  spectrum  the  eigenfunctions  must  then  die  out  exponentially 
at  +00  and  -oo;  therefore,  from  (2.2)  it  is  apparent  that  the  discrete  eigenvalues 
occur  only  for  values  of  k  irtiich  have  non-zero  imaginary  parts.  Since  the  eigen- 

values  of  (2,1)  are  given  by  k  which  is  real  only  irtien  k  is  pure  imaginary  or 

^2 
real  and  since  the  operator  -  — =•  +  V(x)  is  Hermitian,  it  follows  that  the  values 
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of  k  corresponding  to  the  eigenvalues  must  lie  on  the  imaginary  axis.  Since 
when  u(k,x)  satisfies  (2.1)  so  does  u(-k,x),  we  can  without  loss  of  generality 
restrict  the  eigenvalues  to  lie  on  the  positive  imaginary  axis. 

Further,  the  eigenvalues  occur  for  values  of  k  for  >rtiich  the  function 
b(k)  has  poles  on  the  positive  imaginary  axis.  It  is  known  that  these  poles 
are  simple  and  that  b(k)  is  otherwise  regular  in  the  upper  half-plane  and  on 
the  real  axis,  except  perhaps  at  zero.  If  we  add  the  condition  that  xV(x) 
belong  to  L, ,  then  b(k)  has  only  a  finite  number  of  poles  in  the  upper  half- 
plane;  i.e.,  there  are  only  a  finite  number  of  eigenvalues. 

The  completeness  of  the  spectral  resolution  of  the  operator  -  — »•  +  V(x) 


of  (2.1)  can  be  expressed  succinctly  by  the  relation 

00 

(2.9)  2I«(V^  u*(k^,y)  ♦  /  u(k,x)  u*(k,y)  dk  -  5(x-y)  , 

-00 

where  5(x-y)  is  the  Dirac  delta-function  irtiich  has  the  property  that 

00 


S? 


L 


f (x)  5(x)  dx  -  f (0) 

"-00 


for  arbitrary  f (x).  The  function  u  (k,x)  is  a  particular  solution  of  (2.1) 
corresponding  to  the  given  solution  u(k,x)|  the  sum  in  (2.9)  is  taken  over  the 
values  of  k  corresponding  to  the  eigenvalues.  When  u(k,x)  and  u  (k,x)  are  ana- 
lytic functions  of  k  the  relation  on  the  left  side  of  (2.9)  can  be  written  as  a 
single  integral  in  the  complex  plane, 


/  u(k,x)  u  (k,y)  dk  , 


in  ^ich  the  contour  of  integration,  C,  is  a  large  semi-circle  in  the  upper 
half -plane.  The  relationship  between  the  proper  eigenf unctions  ar*l  the  improper 
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eigenfunctions  [u(k,x),  k  in  the  continuous  spectriunj  is  made  clear  by  this  con- 
tour integral  representation. 

3,  Formulation  of  the  main  problem  and  derivation  of.  the  fundamental  integral 
equation 

The  problem  we  are  going  to  consider  is:  given  the  ref ]e  ction  coefficient 
b(k)  as  a  function  of  k  for  an  equation  of  the  form  (2»l)  with  an  unknown  real 
V(x);  to  find  V(x)  as  a  function  of  x.  The  method  we  shall  use  for  this  investi- 
gation is  due  to  Gelfand  and  Levitan'-  J  and  has  been  summarized  by  M,  Levinson'-  -^ 
in  English. 

The  quantity  V(x)  in  (2.1)  can  be  looked  upon  as  a  perturbation  operator 

^2 
which,  when  added  to  the  operator  -  — ^    of  the  equation 

dx"^ 

h2  o 

(3.1)  A-  u^(k,x)  +  k'^  u^(k,x)  -  0  , 

dx 

^2 
produces  the  operator  -  — «•  +  V(x)  of  (2.1).  We  shall  refer  to  functions 

dx 
u  (k,x)  satisfying  (3«l)  as  solutions  of  the  unperturbed  equation, in  contrast 

to  the  functions  u(k,x)  which  are  solutions  of  the  perturbed  equation  (2.1). 
From  this  point  on,  for  simplicity,  we  shall  assume  that  the  function  V(x) 
vanishes  for  x  <  -  •j, where  6  is  some  real  positive  constant. 

The  first  step  in  the  Gelf and-Levitan  method  is  to  recognize  a  funda- 
mental property  of  the  relationship  between  a  solution  u  (k,x)  of  the  unper- 
turbed equation  and  a  solution  u(k,x)  of  the  perturbed  equation.  There  exists  a 
linear  transformation  T(x,y)  which  transforms  a  solution  u  (k,x)  of  the  un- 
perturbed equation  into  a  solution  u(k,x)  of  the  perturbed  equation,  namely, 

00 

T(x,y)  -  /  u(k,x)  u*(k,y)  dk 
-oo 
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where  u  (k,y)  is  the  conjugate  function  of  u  (k,y)  in  the  completeness  relation 


00 


u^(k,x)  ti^(k,y)  dk  -  6(x-y). 


-00 


The  fundamental  observation  to  be  made  is  this:  if  u(k,x)  «•  u(k,x)  for 
X  <  -  5  I  we  recall  that  V(x)  ■  0  for  x  <  -  ^|,  then  T(x,y)  is  triangular 
in  the  sense  that  T(x,y)  «  0  for  y  >  x.  It  is  convenient  to  write 
T(x,y)  «  5(x-y)  +  K(x,y).  Then  K(x,y)  is  also  triangular,  that  is, 
K(x,y)  -  0  for  y  >  x.  We  now  state  the  following: 

Theorem   If  u(k,x)  ■  u  (k,x)  for  x  <  -  r-,then 

x 

(3.2)     u(k,x)  -  u^(k,x)  +  /  K(x,y)  u^(k,y)  dy  , 

-(x+5) 

and  K(x,y)  is  uniformly  bounded,  continuous,  and  has  piecewise  con- 
tinuous first  and  second  derivatives, 

.2 
To  prove  the  theorem  we  apply  the  operator  — »■  -  V(x)  to  (3,2)  on 

dx'^ 
both  sides.  Using  (2.1)  and  (3,1)  we  obtain 


0                          0                                                                     AV(        \                        ^    (k,  x) 
-k'^u(k,x)  -  -kS^(k,x)-V(x)u^(k,x)  ♦  u^(k,x)  ^^^^  ♦  K(x,x)  —2— 

i\y:(-r   -Y  R'*           du  (k,-x-6) 
+  u^(k,-x-5)  ^^^y  ^^  +  K(x,-x-6)  -^—^ *  K^(x,x)u^(k,x) 

X 

+  K^(x,-x-5)u^(k,-x-6)  ♦  /   K^(x,y)u^(k,y)dy 

-(x+5) 

X 

-  V(x)  /   K(x,y)  u^(k,y)  dy, 
-(x+6) 

But  from  (3.2)  and  (3»l),  integrating  by  parts,  we  have 
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.ACkju)  -  -A^(k,x)  +  /   K(x,y)  -^  u^(k,y)  dy 


du  (k,x) 
-k'^u^Ckjy)  +  K(x,x)  -^ ♦K(x,-x-6) 


du  (k,-x-6) 
dx 


•K  (x,x)u^(k,x)  +K  (x,-x-5)  u^(k,-x-6) 


X 


Kyy(x,y)  u^(k,y)  dy 


Then  combining  these  two  results  and  combining  terms  appropriately  we  obtain 


^(x*6) 


K^(x,y)-Kyy(x,y)-V(x)K(x,y) 


Uo(k,y)  dy 


*  h   ^^^  -  V(x) 


-  u  (k.x)  + 
o  * 


dK(x,-x-5) 


u^(k,x)  • 


This  equation  will  be  satisfied  if 


K^(x,y)  -  Kyy(x,y)  -  V(x)K(x,y)  -  0 


and 

(3.3) 

and  if. 


dMx^  -  h  V(x), 


dx     1 


dK(x,-x-5)  _  ^ 


dx 
that  is,  if  K(x,-x-6)  -  constant. 

On  the  other  hand  we  know  from  the  theory  of  hyperbolic  partial  dif- 
ferential equations  that  there  exists  a  K(x,y)  satisfying  (3.3)  (this  is  a 
characteristic  boundary- value  problem).  It  is  easy  to  see  that  having  deter- 
mined K(x,y)  from  (3.3)  we  could  retrace  our  steps  and  arrive  at  (3,2 )»  Thus 


-li- 


the existence  of  K(x,y)  is  proved  (obviously  this  function  satisfies  the  con- 
ditions of  the  theorem).  It  should  be  observed  that,  according  to  (3»3)> 
K(x,y)  is  independent  of  the  choice  of  particular  solutions  u(k,x)  of  (2.1) 
and  u  (k,x)  of  (3.1)  except  that  (3.2)  implies:  u(k,x)  ■  u  (kx)  when  ^  <  ^  -^  * 
Thus  K(k,x)  transforms  any  solution  of  (3.1)  into  a  solution  of  (2.1)  which 
is  the  same  function  for  x  <  -  ^  . 

By  a  completely  analogous  argument  we  can  demonstrate  the  existence  of 
a  triangiJ.ar  operator  kernel  K  (x,y)  such  that 


O.h) 


u^(k,x)  -  u(k,x)  +  /   K^(x,y)  u(k,y)  dy, 
^(x+5) 


when  u(k,x)  •  u  (k,x)  for  x  <  -  »■  • 

In  order  to  make  use  of  contour  integration  in  the  complex  plane 
we  shall  work  with  two  solutions  v(k,x)  and  v(k,x)  of  (2.1)  which  are  ana- 
l3rtic  functions  of  k  for  Im  K  >  0  and  continuous  functions  of  k  for  Im  k  >  0. 
They  can  be  defined  by  their  asjnnptotic  behavior: 


v(k,x)  - 


,a(k)e^^  ♦  p(k)e"^^. 


ikx 


X  •♦  -00 


+  00, 


(3.5) 


A 

v(k,x)  - 


-ikx 


.-plkj  e^*™  ♦  a(k)  e'^^,  x  -»  +oo 


From  the  spectral  representation  theoiry  of  the  operator 
we  have  the  relation 


^ 


*   V(x) 


(3.6) 


k  /   '-i^i^  d.  .  6(x-y)  . 
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where  the  contour  C  is  a  large  seHd-circle  in  the  upper  half  of  the  complex  k- 
plane  and  the  integration  direction  is  clockwise.  The  completeness  relation 

(3.6)  could  also  be  written 

00        ^ 

(3.7)  \^f     ^^^'^ilffiy^  dk-^  a^v(k^,x)v(k^,y)  -  5(x.y)  , 

*-oo 

where  the  sum  rxrns  over  the  poles  of  tttt  in  the  upper  half- plane,   and  the  coef- 
ficients a     are  i  times  the  residues  — ; of       i.  i     at  these  poles. 

Eefore  we  can  proceed  we  must  establish  the  following  identity: 


00 


X  >  y  . 


»-8)  hf    "-^    »'"■"  ^  -  E    V'-v'^''  '■"^'  -  °' 

-00 

To  this  end  we  consider  the  kernel  K(x,y)  in  the  relation 


(3.9)  e"^^-  v(k,y)  +  /   K(y,z)v(k,z)  dz  . 

^(y+6) 

Assume  x  >  y.  Multiply  (3.9)  by  ^^^^^^\  and  integrate  with  respect  to  k 
from  -00  to  +ooj  then  add  to  this  the  result  of  multiplying  (3.9)  with  k  ■  k^ 
by  a  u(k  ,x)  and  summing  over  v.  The  left  side  of  the  result  is  identical 
with  the  left  side  of  (3.8).  By  (3.7)  the  right  side  of  the  result  is 

y 


A 


5(x-y)  ♦  /    K(y,z)6(x-z)  dz  . 
'^(y+6) 

Thus  for  X  >  y  the  right  side  of  the  result  vanishes, and  (3.8)  follows. 

We  are  now  prepared  to  obtain  an  integral  equation  for  K(x,y),  from 

which  by  (3.3)  we  can  determine  V(x).  Consider  the  relation  (3.2)  in  the 

particular  form 

X 

(3.10)     v(k,x)  -  a(k)  e^^  +  p(k)  e"^^  *  /   K(x,z)ra(k)  e^^  p(k)  e"^^^l  dz. 

'^(x*5)      '' 
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--iky 
We  multiply  (3.10)  by    = — HTFT    ^^^  integrate  over  the  contoiir  C  of  (3«6). 

From  (3.8)  it  is  apparent  that  the  left  side  of  the  resulting  equation  van- 
ishes for  X  >  y.     We  have  then  an  equation  of  the  form 

X 

(3.11)  0  -  B(x+y)  +  K(x,y)  ♦  /       K(x,z)  B(y+z)  dz  ,         x  >  y 


^'-(x+e) 


where 


00 


(3.12)  B(r )  -  1^    /     b(k)  e'^^'^dk  -  i  ^  b 

-00 


-ik  r    . 


In  (3.12)  we  have  written  b  for  the  residues  of  b(k)  at  the  poles  k  of  b(k) 

1  /^     ik  f 
in  the  upper  half- plane,  and  we  have  used  the  fact  that  ^  /   e    dk  '  bit) 

-00 

as  well  as  the  definition  b(k)  -  |i^  to  arrive  at  (3.11). 

For  every  fixed  value  of  x  equation  (3.11)  is  a  Fredholm  integral  equa- 
tion of  the  second  kind  for  K(x,y)  when  y  ranges  over  the  interval  -oo  <  y  <  x. 
The  kernel  B(y+z)  is  given  by  (3.12)  in  terms  of  the  reflection  coefficient 
b(k). 

If  (3.11)  has  a  unique  solution  and  we  solve  for  K(x,y),  the  solution 
must  be  the  function  K(x,y)  in  (3.10) •  The  necessary  step  remaining,  therefore, 
is  to  prove  the  uniqueness  of  the  solution  of  (3.11).  We  can  do  this  by  show- 
ing that  the  operator  kernel  6(x-y)  ♦  B(x+y)  is  positive  definite. 

First  we  establish  the  identity 

(3.13)      5(x-y)  .  ^  /  b(k)  e-^^^^-^^dk  -  ^^  /  |a(k)|2  eik(x.y)  ^ 


-00 

o 


*5ir 


^  /  p-  .  b(k)  e-i-] 


.^^  *  b(k)  e-^^ 


dk  . 


-00 

Here  a(k)  is  the  transmission  coefficient*  Using  (2»8)  we  can  write  (3.13)  in 
the  form 


- 11,  - 


-00  o  o 

o  o 

+ 

"-00  -00 

o  o 


1./  ."(^-y'dk  *  jL/  |b(k)|2  e-^^'^-y'dk 


.  1  /  b(ic)  .-^"(^y'  .  Jr  /  bUn  e^''(='*^' 


*/ 


'-00  ''-00 


dk 


The  sum  of  the  first  and  third  terms  on  the  right  side  is  6(x-y),  and  the 
second  and  fourth  terms  on  the  right  side  cancel*  We  have  finally 

(3.1U)  k    ^^^^^   e-i^(-*y>dk  -  ^   /b(k)  e-i^^-^y^dk  .  ^/^m  ei^(^*y>dk 
^00  '^oo  ^00 

Changing  variables  in  (2.U)  and  (2.7),  from  k  to  -k  in  the  second  terms  on  the 

right,  we  see  that  (2.1U)  is,  indeed^ an  identity.  Since  our  steps  are  obviously 

reversible,  it  follows  that  (3.13)  is  an  identity. 

The  term  -i  5~  b  e"^^^^*^  can  be  shown  to  be  real  and  positive  for 
*—  V 


V 

V 


all  real  x  and  y»  The  exponential  factor  is  obviously  real  and  positive,  since 

the  k^  are  pure  imaginary  quantities.  Now  b  -  --,-■-  ,by  the  definition  of  b 

a  (k^) 

as  a  residue  of  a  pole  of  b(k).  From  the  known  relations  (see  [9],  Appendix  II) 

00 
a'(k^)  -  -i  /  v(k^,x)v(ky,x)  dx,  v(k^,x)  -  p(k^)v(k^,x)  , 
-co 

and  from  the  fact  that  v(k  ,x)  is  real  for  real  x,  it  is  seen  without  difficulty 

that  -i  b  is  a  real  positive  number. 

Now  if  there  were  two  solutions,  K-(x,y)  and  KpCxjy),  of  (3.11),  their 

difference  K,-  Kp«  W(x,y)  would  be  a  solution  of  the  corresponding  homogeneous 

integral  equation 
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(3.15)  /  [5(y,2)  +  B(y*z)]  W(x,z)  dz  -  0  . 

■00 

According  to  (3»12)  and  (3»13)  this  can  be  written 


Uy   |a(k)|2  e'^iy^)  dk  .  I,y    [e^^  .  bU)  e"^^] 


[e^^=^  ♦  b(k)  e'Hdk  -  i 


Mx+6)L        "b 


^—     V 


-ik  (y+z)' 


-  W(x,z)  dz  *  0. 


On  multiplying  through  by  W(x,y)  and  integrating  with  respect  to  y  from  -(x+6) 
to  X,  and  then  interchanging  integral  signs,  we  obtain 


00 

h/  I 


a(k) 


.A. 

/  W(x,y)e^^  dy 


-(x+6) 


0    X 

dk  +  ij  /  /  W(x,y)  [e^^+  b(k)e-^^]dy 
^00  ^-  (x+5) 


dk 


(3.16) 


-  r  (-i\) 


iV 


/  W(x,y)  e  "  dy 
*^(x+5) 


0  . 


Since  each  term  on  the  left  side  of  (3.16)  is  positive,  each  term  must  be  separately 


zero*  Thus,  for  example 

X 

W(x,y)  e*^^  dy  -  0, 
"(x+6) 


whenever  |a(k)|  /  0.  It  follows  that  V.'(x,y)  ■  0  almost  everywhere  when 
-(x+S)  <  y  <  x. 


k*     Sufficient  conditions  for  b(k)  to  be  a  reflection  coefficient  (the  synthesis 
problem) 

We  have  seen  how  to  obtain  an  integral  equation  for  the  function  K(x,y), 

from  vriiich  we  can  determine  V(x),  if  the  reflection  coefficient  associated  with 
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V(x)  is  given.  In  the  problem  of  synthesis,  however,  we  are  not  only  required 
to  construct  V(x),  but  generally  we  must  first  fix  upon  an  appropriate  function 
b(k)  for  the  reflection  coefficient,  one  that  satisfies  certain  physical  re- 
quirements. The  question  arises:  if  a  function  b(k:)  is  given,  then  by  the 
process  described  in  Section  3  can  we  always  arrive  at  a  function  V(x)  cor- 
responding to  which  b(k)  is  the  reflection  coefficient?  Undoubtedly,  certain 
conditions  must  be  imposed  on  b(k)  before  this  will  be  possible.  We  shall 
require,  accordingly: 

I,  b(k)  -  b(-k)  for  real  k  ; 
II.  b(k)  is  analytic  in  k  with  no  singularities  in  the  upper  half- 
plane  except  possibly  a  finite  number  of  simple  poles  on  the  positive  imagin- 
ary axis,  and  no  singularities  on  the  real  axis  except,  perhaps,  at  zero; 

III.  The  residues  of  b(k)  at  its  poles  on  the  positive  imaginary 
axis  are  pure  imaginary  numbers  with  positive  imaginary  parts> 
IV.  |b(k)|  <  1  for  all  real  k; 
V,  b(k)  »  e"^'^^  o(i)  in  the  upper  half-plane; 
VI.  the  Fourier  transform  B(^)  of  b(k)  exists  and  is  continuous 
in  -00  <  f  <  00  j 

VII,  BCT)  has  piecewise  continuous  first  and  second  derivatives 
in  any  finite  interval. 

These  conditions  are  sufficient  for  b(k)  to  be  the  reflection  coef- 
ficient corresponding  to  some  quantity  V(x)  in  (2.1),  as  we  shall  see.  Condition 
V  guarantees  that  BCC)  -  0  for  T  <  -6,  and  therefore,  by  (3.11),  that  K(x,y)  -  0 
f or  x  +  y  <  -6.  This  implies  that  V(x)  -  2  ^^^^^^^  -  0  Csee  (3.3 )J  f or  x  <  -  |. 
Conditions  VI  and  VII  along  with  (3»ll)  imply  that  K(x,y)  is  continuous  and 
has  piecewise  continuous  first  and  second  partial  derivatives. 
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We  shall  now  prove  that  b(k)  satisfying  conditions  I  to  VII  is  actually 
the  reflection  coefficient  corresponding  to  some  differential  equation  of  the 
form  (2.1)  with  a  V(x)  obtained  from  b(k)  by  solving  (3.11)  for  K(x,y)  and  set- 
ting V(x)  -  2  — j^^"^  •  The  proof  is  divided  into  two  parts.  In  the  first 
part  we  show  that  K(x,y)  satisfies  (3»3).  We  start  with  (3.11),  whicl^  because 
of  V  [observe  that,  since  x  >  y,  -(y*6)  >  -(x+5)]|,  takes  the  form 

X 

(U.l)  0  -  B(x+y)  +  K(x,y)  +  /   K(x,z)  B(y+z)  dz;   x  >  y,  x  >  -  (y+6). 

*^(y*5) 

Differentiating  (U«l)  twice  with  respect  to  x  gives 

0  -  b"  (x+y)  +  K^(x,y)  +  b' (x+y)K(x,x)  ♦  B(x+y)  ^  K(x,x) 

JUL 

(U.2)  X 

+  B(x+y)  K^(x,x)  ♦  /   K^(x,z)  B(y+z)  dz  . 

-(y*6) 

Differentiating  (U.l)  twice  with  respect  to  y  and  integrating  twice  by  parts 
gives,  when  we  use  the  fact  that  B(^)  is  continuous, 

0  =»  b"  (x+y)  +  K  (x,y)  +  B'(x+y)  K(x,x)  -  B(x+y)  K  (x,x) 


(U.3) 

"-(y+S) 


+  /    ^yy^^»^^  B(y+z)  dz. 


Now  we  subtract  (U.3)  from  (U.2)  and  add  to  this  result  the  product  of 
-2  '^^^^^^  with  (U.l).  We  then  have 
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0  - 


K_(x,y)  -  Kyy(x,y)  -2  ^%^  K(x,y) 


XX 


dx 


(U.U) 


-(y+6) 


K  (x,z) 

XX   ' 


V^.^)  -2^%^K(x,z) 


^x 


B(y+z)  dz. 


It  follows  from  the  uniqueness  theorem  already  proved  that 


(U.?) 


K_(x,y)  -  Kyy(x,y)  -2  ^^^    K(x,y)  -  0; 


XX' 


i.e.,  K(x,y)  satisfies  (3.3 )♦  Further,  a  solution  u(k,x)  of  (2.1)  is  given 
by  (3.2)  when  u  (k,x)  satisfies  (3»l).  The  function  u(k,x)  given  by  (3.2)  is 
evidently  continuous,  and  because  K(x,y)  is  continuous  and  partial  derivatives 
of  K(x,y)  are  piecewise  continuous,  u(k,x)  has  a  continuous  derivative. 

For  the  second  part  of  the  proof  we  must  show  that  the  solution  of 

(2.1)  given  by  (3.2)  when  u^(k,x)  -  e^^  +  b(x)  e"^*"'  is  the  one  which  has 

ikx 

the  behavior  u(k,x)  ~  a(k)  e         for  x  -»•  +oo.     It  is  difficult  to  ascertain 

directly  the  behavior  as  x  -»-  +oo  of  a  solution  u(k,x)  of   (2.1)  when  u(k,x) 
is  given  in  the  form   (3.2),   but  its  behavior  as  x  -»•  -oo  is  just  u   (k,x).     We 
shall  therefore  investigate  the  character  of  our  solution  u(k,x),  which  aj>- 
preaches  e       +b(k)   e  as  x  -►  -oo,    by  comparing  it  with  that    (unique)   solution 

u(k,x)  of   (2.1)  iirtiich  has  the  behavior  e~         as  x  -»•  -oo.     We  know:   if  and 
only  if 


(U.6) 


i    /    u(k,x)  u(k,y)  dk  -  5(x-y)  , 


u(k,x)  has  the  behavior 


(U.7)       u(k,x). 


i 
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e    +  b(k;  e     ,       x  •♦  -oo 
a(k)  e^^        ,       X  ->  +00 


for  the  general  spectral  theory  of  (2.1)  shows  that  (U.6)  is  a  necessary  con- 
dition and  it  follows  immediately  on  multiplying  (U.6)  by  the  conjugate  func- 
tion [the  one  which  has  the  behavior  (U.7)  of  u(k,y3  and  integrating  with  re- 
spect to  y  that  (U.6)  is  also  a  sufficient  condition.  We  write  then 

X 

u(k,x)  -  e^^  +  b(k)  e"^^  +  /   K(x,z)  [e^^^  +  b(k)  e'^^^]   dz 

-(x+5) 
(U.8) 

X 

u(k,y)  -  e-^^+  /    K(x,2)  e"^^^  dz  . 
Mx+5) 

If  we  take  the  product  of  these  two  expressions  and  perform  the  integration 
indicated  in  (U.6),  making  use  of  (3.11),  we  find  that  (U.6)  is  satisfied. 
Thus  the  solution  u(k,x)  of  (2.1)  has  the  required  asymptotic  behavior  and 
the  quantity  b(k)  is,  indeed,  the  reflection  coefficient  corresponding  to 


V(x)»2^%^     . 

5»     Example 

Let  b(k)  -  - 

1 

k2.1 

be  the  reflection  coefficient  for  some  as  yet  un- 
known function  V(x)  in  (2.1).  The  sufficient  conditions  I-VII  on  b(k)  (see  5ec.  U)) 
are  easily  verified  for  this  case.  Also  the  function  B(T)  is  easily  calculated 
by  a  straightforward  application  of  the  theory  of  residues: 


00 


(5.1)       B(T)--i-/  ^ dk  +  ieT 


-ik  sinh  -T,  "^  >  0 


-00  >-        »  I-  _ 
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Then  equation  (3.11)  becomes  in  this  case 

X 


(5.2)       0  -  sinh  (x+y)  +  K(x,y)  +  /  K(x,z)  sinh(y+z)  dz,    x>  |y| 

-y 

We  can  obtain  a  differential  equation  for  K(x,y)  by  differentiating  (5.2) 
twice  with  respect  to  y.  A  single  derivative  gives 


X 


(5.3)       0  -  cosh  (x+y)  +  K  (x,y)  +  /  K(x,z)  cosh(y+2)  dz, 
and  a  second  derivative  gives 

X 

(5,U)       0  -  sinh  (x+y)  +  K  (x,y)  +  K(x,-y)  +  /  K(x,z)sinh(y+z)  dz. 

The  integral  terra  of   (5«U)  is  the  same  as  that   of   (5.2).     We  can  therefore 
eliminate  this  term,   obtaining 

(5.5)  V^""'^^  *  ^^^*'^^  '  "^^^^^^  '  °* 

It  is  evident  from  (5,5)  that  K  (x,y)  is  an  odd  function  of  y.  Thus 

(5.6)  K(x,y)  -  -  K(x,-y)  +  a(x)y  +  2b(x)  , 

where  a(x)   and  b(x)  are  arbitrary  functions  of  x  alone.     Then  (5.5)  becomes, 
after  we  substitute  for  K(x,-y)  from  (5.6), 

(5.7)  K     (x,y)  -  2K(x,y)  -  a(x)y  -  b(x). 
The  general  solution  of   (5.7)  is 

(5.8)  K(x,y)  -  iWl^-2bW  ,     AOO     ,^/2y     ,     BOO     ,-v/2y     ^ 


where  A(x)  and  B(x)  are  new  arbitrary  functions  of  x  alone.  After  substituting 
the  value  of  K(x,y)  obtained  from  (5.8)  into  (3.5),  we  see  that  a(x)  ■  0  and 
A(x)  -  B(x).  Thus,  by  (5.8) 


-  21  - 


iS.9) 


K(x,y)  ■  b(x)  +  A(x)  sinh  \/?  y  . 


To  determine  the  functions  b(x)  and  A(x)  we  substitute  ($.9)  into  ($.2), 
Then  we  have  the  result 


(5.10)      A(x) 


sech  v/?  X            tanh  v/?  x 
.  ,    b(x) 


^ 


^ 


From  (5*9)  we  now  have 


f       -    sinh  y2  X  +  sinh  v^  y 
^  cosh  (/?  y 


(5.11)  K(x,y)  -  . 


y  >  -X 


y  <  -X 


By  (3.3)  we  finally  have 


(5.12)  V(x) 


-  k  sech    y/2     X  ,  X  >  0 


X  <  0     . 


For  this  V(x),  the  solution  u(k,x)  of  (2.1)  having  the  asjrmptotic  behavior  given 


by  (2.6)  and  thus  exhibiting  the  reflection  coefficient  b(k) 
written  immediately  through  the  use  of  (3.2)  and  (5.11) t 


k^.  1 


can  be 


ti(k,x)  -  e^^  -   1    e-^^ 

]£+    1 


sech  \/?  X 


(5.13) 


X 

•  /  [sinh  ^  X  +  sinh  ^  y  j  \e^^-  -^  e"^^]  dy 


After  performing  the  indicated  integration  and  simplifying  the  right  side  of 
(5.13)  we  obtain 
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(5.li4)   U(k,x)  -   - 


k  +  ^/l  ik  tanh  y/2   x 


ikx    1    -ikx 
e   -  -K e 

k^+  1 


ikx 


X  >  0 


X  <  0 


Since,  as  x  ->  +oo,  u(k,x)/>^ 


k^-f  v/2  ik   ikx 


k^+  1 


e    ,  the  transmission  coefficient 


is  given  by 

(5.15) 


a(k) 


k^+  v/2  ik 
k^.l 


Obviously  a(k)  and  b(k)  satisfy  the  required  conditions 

a(-k)  -  aTFr,  b(-k)  -  b(H,   |a(k)|^+  |b(k)|^  -  1  for  real  k. 

We  might  note  in  passing  that  if  V(x)  ■  -  U  sech  \/2  x  for  the  entire 

interval  -oo  <  x  <  oo,  it  is  • transparent ' j  that  is,  there  exists  a  solution 

of  (2.1)  (we  can  obtain  this  solution  from  (5.1U))  which  has  the  behavior 

UA'e    asx->-oo  and  the  behavior  u  i^  3l_— i.  e    as  x  -^  +  oo,  so 

k  -  ,/7  i 

that  no  reflected  wa^e  exists  for  an  incident  wave  from  the  left.  The  ampli- 
tude of  the  incident  wave  remains  one  as  the  wave  approaches  +oo,  but  the  phase 
changes  by  an  amount  depending  on  the  frequency. 


6.  Solution  of  a  general  class  of  problems 

The  method  employed  in  Section  5  of  solving  integral  equation  (5*2) 
suggests  that  a  whole  class  of  problems  associated  with  reflection  coefficients 
of  a  certain  form  might  be  solved  explicitly. 

Suppose  we  have  a  reflection  coefficient  b(k)  satisfying  conditions 
I  to  VII  of  Section  U  and  suppose  also  that  b(k)  is  a  meromorphic  function 
with  a  finite  number  of  poles  in  either  half-plane.  The  function  BCf )  then 
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has  the  form 


(6.1) 


B(t)  -  - 


N      -ik  t 


n-1 


0  »       r  <  -6    • 

Accordingly  the  integral  equation  (3»ll)  for  K(x,y)  becomes 


(6.2) 


-ik  (y+z) 


N      -ik„(x+y)  N       n  -ik  (y+z) 

0  -  E  b  e   "    +  K(x,y)  *  IZ    K    /        K(3f,z)  e   '    dz  . 
n-1  ^  n-1    <//_^c\ 


-(y*5) 


We  can  obtain  a  new  equation  from  (6.2)  by  differentiating  with  respect  to  y| 
each  time  we  differentiate  we  obtain  another  equation.  Let  us  differentiate 
N  times,  forming  N  +  1  equations  altogether,  including  (6.2)  itself.  The  new 
equations  have  the  form 

N  -ik  (x-ty)        /   X  N  ^,       -ik  y    /  -ik  z 

0  »  r  (-ik  )  ^  e      °  ♦  K(^^(x,y)  *  ^  ("i^jV      ^    /    K(x,z)e      "  dz 


(6.3) 


ir«l 


N  ik  5  v-1 


-(r-6) 


lib     e     ""    IZ     ^'^K^^  K^''"^'^^x,-(y*5));         v  -  1,    ...,   N  , 


wl 


.(V), 


th 


where  K^      (x,y)  is  written  for  the  v       y-derivative  of  K(x,y)  and  where  we 

N  ik  5 

have  used  the  fact  in  (6.1)  that  B(-5)  -  T~    be     "    "0.     We  can  eliminate  the 


n-1 


N  quantities 


/   K(x,z)  e   "  dz 


>,th 


from  equations  (6.3)  and  (6.2)  obtaining  a  single  N  -order  differential  equa- 
tion in  K(x,y)  and  K(x, - (y+6 ))  of  the  formt 

N-2 


(6.U)   0  -  P(x,y)  +  21  V^^  K^  ^x,y)  *  ^1    V^^  ^    (x,-(y+6))  . 
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By  substituting  -y-5  for  y  in  (6«U)  we  obtain  another  differential  equation, 
3y  differentiating  (6.1i)  and  then  performing  this  substitution  again  we  obtain 
a  new  equation.  If  we  differentiate  N  -  2  times  and  perform  the  substition  each 
time,  we  arrive  at  a  system  of  equations  involving  K(x,y)  and  its  derivatives 
and  K(x,-(y+5))  and  its  derivatives.  It  is  now  possible  to  eliminate  K(x,-(y+5)) 
and  its  derivatives  and  finally  obtain  a  single  differential  equation  of  order 
2N  -  2  for  K(x,y).  In  the  example  in  Section  5,  N»2,  and  consequently  we  obtain 
a  second-order  differential  equation  for  K(x,y)  in  this  case. 

In  most  cases  we  can  solve  the  integral  equation  (3»ll)  by  a  more  direct 
method,  however.  Let  us  assume  that  K(x,y)  has  the  form: 

N         ay 

(6.5)  K(z,y)  -  n   ^v^^^  «     » 

V— N 

where  the  sum  extends  over  all  integers  between  -N  and  N  except  zero,   and  where 
a     and  f   (x)  are  quantities  to  be  determined.     If  we  substitute   (6.5)  into   (3.11) 
we  have 

N  -ik  (x+y)      N  ay       N       N  -iky  p     (a  -ik„)2 

0-     r    b^e       °  .n      f^(x)e^     ^^H      b„e       "f,(x)/e     ^       "    dz 


N  -ik  (r»y)      N  a  y     N      N  -ik  y 

n-l  V— N  nil    V— N 


-(y+5) 
(o^-ik^x     .(y^)(a^-lk^) 


-e 


a  -ik 
V       n 


a  X 

N  -ik„(x*y)      N  -ik„(x+y)    N         e  ^  f^(x)         N  ay 

H=l    "  fe  "  t^-N        \^       v=:n    ^ 

ik  6 

N          -a  y-a  6                N        b     e     " 
T—  V      V      -    /    N  r—         n 

v«-N  n-l       V  n 
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This  form  suggests  that  we  assume  a  =  -a  ;  then 

V        -v' 


(6.6) 


N 


0 '  1 +  Y2 


a  X 

v»-N        V  n 


n  «  Ij  • .  • ,   Nj 


and 


(6.7) 


a  6      N 
0  =  f,.(x)  +  f_„(x)  e  ^       J2 

n"l 


-V 


b     e 
n 


ik  5 
n 


a  +  ik    ' 
V         n 


1,...,   N 


A  similar  set  of  equations  goes  with  (6.7)  and  can  be  deduced  from  it  by 
substituting  -v  for  v.     In  order  to  get  a  non-trivial  solution  f   (x)  of  the 
resulting  homogeneous  system, the  coefficient  determinant  must  be  set  equal 
to  zero: 


(6.8) 


ik  6 

N        be" 

rn 

n-1        n        V 


ik  5 
n 


N  be" 

c-        n 

^  ik     +  o 

n=l  n        V 


-  1  -  0  . 


Equation  (6.8)  reduces  to  an  N  -degree  algebraic  equation  for  at   thus  in 

general  there  are  N  values  which  a  can  have.  It  is  evident  from  (6.8)  that  if 

a  is  a  solution  then  a  ■  -a  is  also  a  solution.  When  a  _  is  known,  equations  (6.7 ) 
V  -V     V  V  ' 

can  be  used  to  give  f  (x)  in  terms  of  f  (x).  After  substituting  this  result 
in  (6.6)  we  have 


(6.9) 


N 

-i-H 

v«l 


^  /  ^c\  ik  5 

-^V^^*^^   N    be  ° 

En 
-   IF^ 


a  X 

a  -  ik    a  +  ik    — ^ 
V    n    V     n   n^l 


a 
n   V 


f^(x)|  n  «  1,.,.,  N. 


Equations  (6.9)  can  be  solved  for  fy(x)  provided  that  the  coefficient  determin- 
ant does  not  vanish.  The  coefficient  determinant  will  surely  vanish  if  any  two 
values  of  a  are  the  same  as  v  goes  from  one  to  N.  Thus  we  need  the  condition 
generally  that  (6.8)  have  no  multiple  roots,  except  perhaps  a  »  0^ which  can  be 

a  double  root. 
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If  it  turns  out  that  equations   (6.9)  have  no  solutions,  we  must 
assume  a  different  forin  for  K(x,y): 

N  ay 

(6.10)  K(x,y)  -     ^        fix)  P  (y)  e  "     , 

v»-N      ^ 

where  the  P  (y)  are  polynomials  in  y.  A  similar  procedure  to  the  one  just 

described  should  lead  to  a  determination  of  the  a  ,  f  (x)  and  the  P  (y)  and 

thus  of  K(x,y)» 


7,  The  second  differential  equation 

Often  the  differential  equation  of  interest  in  scattering  problems  is  not 
(2«l)  but  rather  an  equation  of  the  form 

(7.1)  ^^^    *    kV(x)  u(k,x)  -  0  , 

where  n(x)  is  real,  continuous  and  has  a  continuous  derivative  in  -oo  <  x  <  oo. 
The  solution  of  (7»l)  can  also  be  interpreted  as  a  wave  if  n(x)  approaches  a 
constant  as  x  -»  +  oo,  and  a  reflection  coefficient  b(k)  can  again  be  defined. 
The  probLera  of  determining  n  (x)  from  the  reflection  coefficient  can  be  dealt 
with  by  a  consideration  of  an  equation  of  the  form  (2.1)  corresponding  in  a 
certain  way  to  (7.1)« 

Suppose  that  we  have  a  differential  equation  of  the  form  (2.1): 

(7.2)  w"  (z)  +  [k^  -  V(z)]  w(z)  =  0  , 

where  the  primes  indicate  differentiation  with  respect  to  z.  Let  us  make  the 
substitution  x  «  x(z)  in  (7*2).  We  have 

(7.3)  w  ♦  ^  w  +  iiLzVl  ^  .  0, 

X  X 

where  again  the  primes  indicate  differentiation  with  respect  to  z  and  the  dots 
with  respect  to  x.  Now  let  us  make  the  substitution  w  ■  ^^'^^-'^■^  in  (7 •3): 
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(7.h) 


U  + 


3x 


n  2 


4^  ^ 

X      Ux 


2x 


V 
•2 


u  ■  0  . 


If  the  second  term  in  the  bracketed  factor  of  u  in  (7.U)  is  zero,  (7.l4)  will 
have  the  form  (7«l)  with  n(x)  ■  — j-  .  Moreover,  if  x(z)  i:|  a  monotonic  function 

X 

of  z  and  independent  of  k,  the  reflection  coefficient  is  the  same  for  equations 
(7,2)  and  (7«li)«  Therefore,  we  can  find  n(x)  -  ■=j-  if  we  can  find  x(z).  We 

X 

assume  that  V(z)  is  constructed  by  the  method  of  the  previous  sections,  and 
now  it  can  be  used  to  obtain  x(z). 

We  set  the  second  term  in  the  brackets  in  (7«U)  equal  to  zero,   and 
this  implies 

"2  •      "'       .     •2      /    X        -. 

3x  '^     -     2x     X       -  Ux       V(z)  -  0 


or 


or 


(7.5) 


t  in  ti 

XXX 


2^(logx")    -3^(logx')-    3 


-  U  V 


dl  ^^°S  ^   ^ 


tf 
x         d 


Now  we  define  co  =  — j-  =  3~"  (log  ^  )>   from  which  we  have 
(7.6) 


dl  (log  X     )  -    <o  +  --     . 


Substituting  in  (7.5)  gives  us 


(7.7) 


»        1     2 
CO     -  i  CO     -  -2  V  . 


Equation  (7.7)  is  a  Ricatti  equation  and  can  be  written  as  a  linear  second - 
order  differential  equation  by  means  of  the  transformations 
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We  have  finally 


n 


(7.8)  U  (z)  -  V(z)  U(z)  »  0, 


where  x  =  -s and  C  is  an  arbitrary  contant;  thus 

ir(z) 


(7.9)  X  -  C 


z 

dz 


?(z) 


After  we  have  found  V(z)  by  the  methods  of  the  previous  sections,  we 
can  use  (3.2)  with  k  set  equal  to  zero  to  obtain  a  solution  U(z)  of  (7«3). 
Then  by  means  of  (7.3)  and  (7.9)  we  get 

(7.10)  n(x)  -  -T-i r-  -  U^(z(x))  . 

X  (z(x)^ 

The  case  where  n(x)  *  constant  in  (7.1)  corresponds  to  the  case  where 
V(z)  ■  0  in  (7.2).  According  to  (7.9),  x  would  then  be  proportional  to  z  or 
to  -rr»  where  A  is  an  arbitrary  constant.  The  second  case,  however,  leads  in 
general  to  a  singular  n(z)» 

In  each  problem  it  must  be  decided  which  of  the  solutions  of  (7.8)  to 
use,  according  to  the  physical  situation.  Obviously,  n(x)  is  not  determined 
uniquely  by  the  reflection  coefficient  as  is  V(z). 
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